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The one dimensional spin system consisted of triangular 5=1/2 XXZ Heisenberg 
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state phase diagrams, corresponding to different regions of the values of system 
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I. INTRODUCTION 

The Heisenberg model, on the one hand, is the key component of microscopic physics underlying 
the magnetic properties of materials on the other hand, its one-dimensional variant plays one of 
the fundamental roles in mathematical physics as the example of integrable quantum many particle 
systems {2!. However, applicability of its exact solutions (which exists only in one- dimensional 
case) in describing the thermodynamic properties of real magnetic materials are rather limited, 
because of fact that in given context the term "exact solution" means rather possibility to obtain 
quantum-mechanical spectrum of the system in principle, than the possibility of obtaining analytic 
expressions for thermodynamic functions. So, the Bethe ansatz technique, which is applicable in 
simplest case of Heisenberg spin chains, leads to the transcendent equation for determining the 
spectral parameters [3], but it is almost useless in calculating the partition function. Though, there 
exists rather sophisticated analytic technique for calculating the finite-T properties of Heisenberg spin 



chains [4], it is still of current importance to find alternative relatively simple approximate methods 
of describing the thermodjTiamic properties of the model of real magnetic materials. One of such 
methods consists in the changing of some (or even all) exchange Heisenberg interactions between 
the spin with more simple "classical" ones, which allows one to develop the exact classical transfer- 
matrix technique for obtaining analytic expression for partition function and all thermodynamic 
functions. This approach has been shown to be in satisfactory qualitative(even quantitative in some 
cases) agreement with the experimental data and numerical calculations for some variant of one- 
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FIG. 1: The schematic picture of the Ising-Heisenberg chain with triangular plaquettes alternating with 
single sites. The filled(empty) circles indicate Heisenberg(Ising) spins. 



dimensional spin chains, namely in F-F-A and F-F-A-A alternating chains |5|]-[l0|. In this work we 
consider the Ising-Heisenberg chain with triangular Heisenberg XXZ clusters alternating with single 
Ising spins, the interaction between spins in each triangle and its adjoint two Ising spin is assumed 
to be of Ising type which allow one to implement the exact analytic calculation of all thermodynamic 
properties. 1 

II. SOLUTION OF THE GENERAL MODEL OF HEISENBERG-ISING CHAIN 

Let us consider the system of Heisenberg and Ising spins residing at the one-dimensional chain 
in the following way. Heisenberg spins form clusters with relatively small number of sites, within 
these cluster the spins are interact with exchange interaction, another kinds of interaction, e.g. 
Dzyaloshinskii-Moriya term, single-ion anisotropy, dipole-dipole interaction, Zeeman terms, e.t.c. 
can be presented as well. The topology of links between the Heisenberg spins can be various, in 
most simple cases it could be finite linear chain or close loop. The Heisenberg spin clusters compose 
the chain by alternating with the single sites with Ising spins. This means that quantum spins from 
different cluster do not interact to each other, but only with the Ising spins between them. Interaction 
between Heisenberg and Ising spins is supposed to be of Ising type, i.e. S^a. The Hamiltonian for 
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such a model can be written in the following way 

N 



Ti- = ^ yii (Sji, ^im) — 2^ (c^i + j ) 
i=l ^ ^ 



(1) 



where Tii (Sji, Sj^) is the Hamiltonian of the i — th Heisenberg cluster, containing m sites, (Tj and 
(Tj+i are two its adjacent Ising spins, takin value ±1 and H2 stand for the dimensionless magnetic field 
acting on Ising spins. One can assume for a while that the corresponding field acting on Heisenberg 

spins via —HiS^^ differs from the Hi which physically can be explained by the difference in the 
gf-factors. The term describing interaction between the Heisenberg spins of i-th cluster and Uj and 
(Tj+i is included into Hi and in general case has the form: 

m 

'^i" = E ('^^^^ + ^la. (2) 

a=l 

where La{Ra), a — 1, ...,m is the set of the coupling constants of the interaction between a-th spin 
from the cluster and its left(right) adjacent Ising spin. The absence of the immediate exchange 
interaction between spins from different clusters leads to commutativity of the Hamiltonians for 
different clusters 

[ni,nj]^o. (3) 

This fact makes possible to use the classical transfer-matrix technique for obtaining the exact ther- 
modynamic solution for such kind of spin chains. Namely, consider the partition function 

Z = ^TY(5)e-^^-i(«'-'^('^'+^'+^)), (4) 

where the sum is going over all possible configurations of the Ising spins and Tr(5) denotes the trace 
over all Heisenberg operators 5" and /3 as usually is the inverse temperature. Commutativity of the 
Hamiltonians for different cluster allows one to expand the exponent in the partition function: 

N N 

(a) i=l (a) i=l 

D D 

n=l n=l 

here TVj stands for the trace over the state of i-th Heisenberg spin cluster and D — YYk=i{'^^k + 1) is 
the dimensionality of the corresponding matrix, sj. is the spin of k-th spin in the cluster. Thus, the 
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partition function takes the form of the partition function of one-dimensional chain with classical 
discrete two-state variable on each site. The corresponding transfer-matrix T reads : 



(6) 



where Z+ = Z{1, 1), Z_ = Z{-1, -1), Zq = Z{1, -1) and Zq = Z{-1, 1). According to Eq. (El), the 
partition function is the sum of N-th powers of the eigenvalues of T, which are given by: 

= 1 (^e^^^Z+ + e-^^2Z_ ± ^ {e^H^Z+ - e-P^^Z^f + AZoZ^ . (7) 

Thus, the free energy per one spin in the thermodynamic limit when only maximal eigenvalue survive 
is 

/ = - ^^l^^^ log Q (e^^^ + e-^^^ Z_ + sJ{eP^^Z^-e-P^^Z.f + AZoZ^ ) , (8) 

The different variant of Ising-Heisenberg chains with two- and three-spin Heisenber g cl usters with 
the properties described above are considered in Refs. 

Heisenberg chain with triangular XXZ spin-1/2 clusters was solved exactly within the technique 
described above. The system is described by Hamiltonian given by Eq.([l]) with 

+K{St, + St2 + St,){a, + - //^l^A + + St,), 

All free spins from each triangle interact with their two adjacent Ising spin with the same constant 
K (See Fig. ([T])). To construct the transfer matrix T one should obtain eight eigenvalues An(o"j, cxj+i) 
of these Hamiltonian, which are the following 

Ai,2(fx„ = 3(J ±H)t 3K{(Ti + (Xi+i), (10) 

A3,4(cri, CTi+l) = A5,6(cri, CTi+i) = -(1 + 2A) J ±HT K{(Ti + CTi+i), 



Thus 



Z{(Ji, a.+i) = 5i cosh(/5(i/2 - K{ai + a.+i))) + B2 cosh(/33(if2 - K{a, + (11) 



Having all these functions one can easily obtain the partition function and free energy and then, 
exploiting general thermodynamic relations, taking a derivatives of it, various thermodynamic func- 
tions, like magnetization, specific head, susceptibility, e.t.c. It Ref. ([ly]) these issue was presented 
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in details for the case of antiferromagnetic coupling. Here we analyze the ground state properties of 
the system in all range of coupling constants. 



III. T=0 PHASE DIAGRAMS 

Let us write down the possible microscopic configurations of spins in the chain with corresponding 
energies per one lattice unit (in our case each block contains one triangle and its one adjoint Ising 
spin) and magnetizations. Up to the inversion of the direction of total magnetization, one can obtain 
six different eigenstates of the model under consideration: saturated completely polarized state (S), 
three ferrimagnetic states, with magnetization M = 1/2, (-Fi), (-^2); (-P3) and two antiferromagnetic 
states with M = 0, (AFi) and {AF2): 

N 

|5) = ni3/2,3/2),(g)| T)„ Es = 3J + 6K-4H, M = 1, (12) 



i=l 

N 



Fi) = n 13/2,3/2) 
1=1 

N 

F2) = n 11/2,1/2) 
1=1 

N 

F3) = n 13/2,1/2) 

i=l 
N 

AFi) = n 11/2,1/2) 

i=l 
N 

AF2) = n 13/2,1/2) 



Ef, =3J-6K- 2H, M = 1/2, 
t)i, Ef, = -J{1 + 2A) + 2K-2H, M = 1/2, 
t)i, Ep, = -J(l - 4A) + 2K-2H, M = 1/2, 



i). 



Eaf, = -J{1 + 2A) -2K, M = 0, 
- J(l - 4A) -2K, M = 0, 



E 



AF2 



1=1 



where |/,m) stand for the standard symmetry adopted states of three S=l/2 spins with 5*°* = / 
and 5**°* = m and | t (i)) are up (down) states of the corresponding Ising spins. Compa ring the 
energies of different states one can obtain T = phase diagram of the model. In Ref. lOj] the 
diagram was found only for antiferromagnetic values of all coupling constant, here we present the 
phase diagrams for all possible values of model parameters. In Fig. ([21) one can see the ground states 
phase diagram in the {H — A)-plane for the case J > 0, K > and particular value r] = \ J / K\ = 1. 
In this case all six ground states can be realized depending of the value of exchange anisotropy 
A and external magnetic field. For positive values of A one can see (S"), (-Fi), {F2) and {AFi) 
states, whereas for A < (-F3) and {AF2) states appear instead of {F2) and {AFi) respectively. The 
boundary between {AFi) and (Fi) states is given hj H = J(2 + A) — 2K, between {AF2) and (Fi) 
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hj H = 2J(1 + |A|) - 2K, between (S) and (Fa) hj H = J(2 + A) + 2K and finally between {S) 
and (F3) hj H = 2J(1 + |A|) + 2K. The ferrimagnetic phase Fi passes to the saturated phase at 
H = 6K, another horizontal line H = 6K in the diagram separates antiferromagnetic phases from 
ferrimagnetic ones. In case of J < 0, i^' > (Fig. Q) again one can see all six phases. The qualitative 
difference from previous case is the appearance of the wide gap between (AFi) and {AF2) states. It 
is also should be noted that the phases which were compatible with positive values of A in previous 
case now correspond to its negative values and vice versa. Like in the previous case two horizontal 
lines a.t H = 6K and H = 2K separate fully polarized phase from (Fi) and antiferromagnetic phases 
from ferrimagnetic ones respectively. The boundaries between phases are defined by the following 
equations: H = 2\J\{A-1) - 2K between(Fi) and(AF2); H = \ J\{\A\ - 2)-2K between (Fi) and 
(AFi);|A| = ^^^+1-^1^ between (Fi) and (F2); and A = ^^±^ between (Fi) and (F3). The next case, 
J > 0, < (Fig. (jl])) is characterized by only three phases to appear, namely, (F2), (F3) and {S). 
In this diagram two different ferrimagnetic phases are separated by A = value. So, at negative 
values of A the ground state of the system is (F3), positive values of A lead to the (F2) ground state. 
Both ferrimagnetic phases pass immediately to saturated one at certain value of the magnetic field. 
For (F3) phase these value is given hj H = 2J(1 + |A|) - 2\K\, for (F2) hj H = J{2 + |A|) - 2\K\. 
In Fig. ([5]) the ground state phase diagram for J < 0, -ft' < is presented for 77 = 2. Here again only 
three phases from previous case can be realized. For the values of axial anisotropy A belonging to the 
interval from A = — 2-^^^jp^ to A = -^^^jjp^, the only possible ground state of the system is (S). For 
A < — 2^^^jjp^ the ground state became (F2) which pass to saturated one at H = | J|(|A| —2) —21^1. 
The same picture takes place for A > ^"^^j j^^*^^ with respect to the (F3) phase. 
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FIG. 3: r = ground state phase diagram ior J < 0,K > and r/ = 1. 



